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Abstract 

Experiments of Electron Spin Resonance (ESR) were performed on 
Co 2+ substituting Zn 2+ or Mg 2+ in powder samples of Zn2(OH)P04 and 
Mg2(OH)As0 4 . These two compounds are iso-structural and contain the 
Co 2+ in two environments of approximately octahedral and trigonal bipyra- 
mid structures. The observed resonances are described with a theoretical 
model that considers the departures from the two perfect structures. It is 
shown that the resonance in the penta-coordinated complex is allowed, and 
the crystal fields that would describe the resonance of the Co 2+ in the two 
environments are calculated. The small intensity of the resonance in the 
penta-coordinated complex is explained assuming that this site is much less 
populated than the octahedral one; this assumption was verified by a molec- 
ular calculation of the energies of the two environments, with both Co and 
Zn as central ions in Zn2(OH)P04. 

76.30.Fc,71.70.Ch,31.15.-p,31.15.Ct 
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I. INTRODUCTION 



Mineral solid state chemistry offers an important contribution to materials science [IJ in 
the search for systems with new and useful physical properties. The phosphate and arsenate 
minerals crystallize in various structures, sometimes containing several non-equivalent sites 
for the metals. The minerals of the olivine group, with the ABXO^ formula, have been known 
for a long time ||, and the adamite family, with formula [M 2 (0/OH)(AT)4)], belongs to this 
group and takes its name from the natural compound |5|,[3J Zn 2 (OH)As0 4 . The compounds 
studied in this work belong to this family, and the cations can occupy two sites with rather 
different environments, one being octahedral and the other penta-coordinated, so that rather 
different magnetic properties could be expected when magnetic cations are employed. The 
recently synthesized compounds Zn 2 (OH)P0 4 §, Co 2 (OH)P0 4 §, Mg 2 (OH)As0 4 as 
well as the natural Co 2 (OH)As0 4 |L|, are °f the adamite type and present these two type of 
sites. We have then found interesting to study the properties of the Co 2+ ions as impurities 
in the two non-magnetic compounds, as a first step in the understanding of the properties 
of the concentrated compounds. 

Electron Spin Resonance (ESR) is a fruitful technique to obtain local information of the 
environment of the magnetic ion, and the Co 2+ ion is particularly useful for this kind of 
study because its g value has a strong crystal field dependence in these compounds. To 
analyze the ESR measurements it is necessary to have information about the splitting of 
the energy levels with both the crystal field and the electronic Coulomb repulsion, and to 
obtain this information we employed optical diffuse reflectance measurements ||. 

The experimental ESR powder spectra of Co 2+ impurities in both Zn 2 (OH)P04 and 
Mg 2 (OH)As0 4 present two different sets of lines, one very intense, and the other just ob- 
servable. The average of the g factors of the intense spectra is 4.15 in the two compounds, a 
value close to the 4.33 expected for Co 2+ in moderately distorted octahedral symmetry j9], 
and it seems reasonable to assign these spectra to that environment and apply the approach 
already employed |10[ in the study of the ESR of Co 2+ in NH4NiP04.6H 2 , where the 
crystal fields that reproduce the observed spectra were obtained. As the remaining lines are 
very weak one should analyze whether they belong to the penta-coordinated symmetry, and 
in that case the possible reason for their small intensity. 

We are not aware of any theory describing the ESR of Co 2+ in the penta-coordinated 
environment, a distorted trigonal bipyramid, and we first calculated the crystal fields of the 



perfect trigonal bipyramid following the existing literature [11,12]. To analyze the distorted 



complex we derived the normal modes of the trigonal bipyramid with respect to the reference 
complex, and then obtained the Jahn- Teller contributions |T3"l , |T4l to the crystal field acting 
on the Co 2+ , that is generated by these modes. 

In this calculation we have introduced a procedure that uniquely defines the orientation 
and size of the two reference complexes, so that the normal modes that describe their 
deformation are free from irrelevant rotations and expansions. 

These results were then employed to calculate the theoretical ESR spectra. We found 
that for the system parameters obtained from the optical spectra we should expect that the 
ground doublet be Mj = ±1/2, corresponding to an allowed spectrum. The rather small 
intensity, of this type of spectra seems to indicate a preference of Co 2+ for the octahedral 



sites in the crystal structure, a conjecture that was advanced in a preliminary report [15 
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on the ESR of impurities of this ion in Mg2(OH)As0 4 , and that is confirmed in the present 
work. Employing a molecular calculation we have also verified that the formation energies 
of the two type of complexes, with both Co and Zn as the central ions, are compatible with 
this hypothesis. 

The description of the experiments, and experimental analysis of the data is presented in 
Section ||. The theory of the ESR of the Co 2+ in a distorted trigonal bipyramid is presented 
in Section |TI|, together with the theoretical analysis of the two type of complexes in the 
two compounds. A discussion of our results is presented in Section |V] together with our 
conclusions. 



II. EXPERIMENTAL 



A. Synthesis and characterization of the materials. 

Compounds with Co 2+ substituting Mg and Zn in Zn 2 (OH)P0 4 and Mg 2 (OH)As0 4 were 
prepared by hydrothermal synthesis, starting from the (M, Co)3(X04) : 8H 2 (M = Zn, 
Mg) vivianites, previously prepared as reported elsewhere ]TB|. Approximately 0.200 g. 
of these precursors were disaggregated in ca. 75 mil. of water and were placed in a 
poly(tetrafluoroethylene)-lined stainless steel container (about three quarters full) under 
autogenous pressure. The reaction was carried out at 180° C and maintained for one week. 
The resulting micro crystalline products were filtered off and washed with ether and dried in 
air. 

The results of the analysis of Mg, Zn, Co, P, and As by inductively coupled plasma atomic 
emission spectroscopy (ICP-AES) are in good agreement with the proposed formulae. The 
compounds were also characterized by X-ray powder diffraction, using the Rietveld method. 
The diffractograms were indexed with the P„„ m space group and the lattice parameters 
a = 8.042(3) A, b = 8.369(2) A and c = 5.940(2) A for the phosphate compound and 
a = 8.286(2) A, b = 8.594(2) A and c = 6.051(1) A for the arsenate. The parameters of the 
phosphate compound are only slightly different from those given in reference [|5] , while those 
of the arsenate coincide with the data published by other authors ||[7| ; the results obtained 
in the last three references were performed in single crystals. The X-ray powder pattern was 
recorded employing Cu Ka radiation with a PHILIPS X'PERT automatic diffractometer, 
with steps of 0.02° in 20 and fixed-time counting of 1 s in the 5 < 2G < 70° range. We 
preferred to use our own experimental parameters in the present paper. 



B. Optical studies. 



The necessary optical data were obtained from diffuse reflectance experiments, performed 
in a CARY 2415 UV-VIS-IR spectrometer, controlled with a VARIAN DS15 workstation, in 
the 5000 - 50000 cm -1 wavenumber region flB,I7|. The whole of the optical data used in this 



work was recorded at room temperatures, and all the relevant data that was necessary in the 
present work is in the tables | and |TT]. Figure [1] presents the experimental data for both the 
phosphate and the arsenate Co compounds studied in this paper. The system parameters of 
the octahedral complexes are slightly different from those already published P,[PT| , because 



3 



they were obtained from the optical spectra after including a spin orbit correction in the 
ground orbital level [18]. 



C. Electron Spin Resonance (ESR). 

The ESR spectra were performed at X Band on a Bruker ESP300 spectrometer. Cooling 
and temperature control of the samples were obtained with a standard OXFORD helium 
continuous-flow cryostat, included in the microwave cavity. Magnetic field measurements 
were done simultaneously with the ESR spectra recording, using a Bruker ER035M NMR 
gaussmeter. The resonant frequency of the cavity was measured with a Hewlett-Packard 
5352B microwave frequency counter. 

Only powder spectra could be measured for the two systems studied here because it was 
not possible to obtain single crystals, and small concentrations of Co (1% in the arsenate 
and 0.1% in the phosphate) substitute the metals in the two lattices. The curves denoted 
with (a) in Figs. ||] and |3] show the measured ESR spectra for the two samples, recorded at 
4.2 K, and they both clearly show three sets of lines with a well defined hyperfine structure 
that identifies the Co 2+ ion. There are also some extra lines, rather weak in the phosphate 
but more intense in the arsenate, that preclude an automatic fitting of the spectra. We have 
then simulated the powder spectra of the hexa-coordinated Co 2+ with a program which 
allows any symmetry, line positions, hyperfine tensor, and linewidth anisotropies, and our 
best results, plotted in the curves (c) of Figs. |2| and correspond to the g-values shown in 
rows a) of table |T| Their values and positions are also shown by arrows below the simulated 
curves (c). 

The extra lines near 200 mT in the phosphate show an hyperfine structure typical of the 
Co 2+ , and are given in more detail in the inset of Fig. ^|. The remaining lines in the two 
compounds are rather broad and show a collapsed hyperfine structure. The curves (b) in 
the Figs. || and|| show the sum of the simulated spectra of the hexa-coordinated Co 2+ in (c), 
plus a simulation of the penta-coordinated Co 2+ that employs the g- values given in rows b) 
of table |H and is adequately renormalized to account for the smaller relative concentration 
of the last compound. These g values have rather large errors, and their positions are shown 
by arrows above the measured spectra (a) of Figs. |2| and |3|. In the inset of Fig. ^| it is also 
shown the detail of the hyperfine structure near 200 mT both in the experimental and in 
the simulated spectrum. 

The assignment of the extra lines to the penta-coordinated complex shall be further 



discussed in section III B 5 



III. THEORETICAL DISCUSSION 

A. Hexa-coordinated Co 

We shall first discuss the hexa-coordinated Co 2+ ions, that are surrounded by six oxygens 
in a fairly regular octahedron with positions given in table 0. In the present case we have 
a powder spectra, and we could only measure the three principal values of the g tensor 
(see table PTI] ). As in a previous work [[H]] we shall consider the effect that the crystal field 
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generated by the normal modes of the octahedron has on the gyromagnetic tensor g. This 
method systematizes the procedure, and in table [V| we give the normal modes that reproduce 
the experimental values of the three g^. We shall then choose a reference perfect octahedron 
centered in the Co 2+ , calculate the normal modes corresponding to the crystallographic 
positions of the O, and compare them in table [V] with those obtained from the experimental 
spectra. 

Only the normal coordinates that are invariant against inversion with respect to the 
center of the octahedron are necessary in the present problem, [|19"l , |10f and these are separated 
into the three sets {Qi}, {Q2, Q3} and {Q4, Q5, Qq}, with the corresponding Qj transforming 
respectively like the basis of the irreducible representations A\, E and T 2 of the cubic group, 
as given in table II of reference |20|. 



The 4 F ground state of isolated Co 2+ (3d 7 ) in a purely octahedral crystal field splits into 
two orbital triplets 4 T 1; 4 T 2 and one orbital singlet 4 A 2 . Spin-orbit effects partially lift the 
degeneracy of the 4 Ti triplet into one T 6 , two r 8 and one T 7 subspaces, and the resonance for 
the lowest doublet (T 6 ) is isotropic with g = 4.33 0. The addition of lower symmetry crystal 
fields produce further splitting of the 4 7\ triplet, giving six Kramer's doublets, and in most 
cases it is found that the trace of the g tensor is close to the cubic isotropic value; [^J in the 
present case the average g is 4.1537 for Co:Zn 2 (OH)P0 4 and 4.153 for Co:Mg 2 (OH)As0 4 . 
To understand these values it is sufficient to consider the Co 2+ in pure octahedral symmetry, 
because the crystal fields of lower symmetry do not change this value in our approximation. 
The calculation follows the same lines given in reference [[Uj] and shall not be repeated here. 



In the lowest order one obtains g from the matrix elements of the Zeeman term in the Tq 
subspace of the 4 Xi ground triplet. The matrix elements of the orbital angular momentum 
L within a Ti subspace are proportional to those of a P term, but one should note that the 
excited term 4 P is also of the 4 Ti symmetry, and is mixed by the cubic field with the 4 7\ 
of the ground 4 F term. If we indicate two states of 4 F and 4 P with <pi and respectively, 
such that they transform in the same way under the cubic group, the states of the ground 
4 Ti will be of the form a<pi + bfy . The values of the constants a and h can be obtained 



2!||23H from the Racah parameter B and the crystal field parameter D q , that were estimated 
T8| from the spectroscopic data and are given in table |. With these values one obtains 

a = —0.9820 and b = 0.1886 for the phosphate, and the proportionality constant of the 

angular momentum is then 

a = -1.5 a 2 + b 2 = -1.4110. (1) 

To analyze further the experimental g tensor, one could try and find crystal field values 
that would reproduce the measured results, and a study of this type was presented by 
Abragam and Pryce for the Cobalt Tutton salts. To simplify the study we present a 
model that describes all the crystal fields acting on the Co as originating in the crystal 
field of the six nearest O located at the vertices of a deformed octahedron, obtained by 
displacement of the vertices of the reference octahedron. If one neglects the mixing of other 
configurations into the ground configuration (3d) 7 , it is sufficient to keep only the part of 
the crystal field V that is even against inversion. We could then write V = Ya=i ^( r t) > 
where V(r) would be the sum of homogeneous polynomials of second and fourth order in 
the components x, y, z of the electronic coordinates r. Within our model, one could then 
write [1191 
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V(r)=£Q^.(r) 

3 



(2) 



where the Qj and Vj(r) transform like the same partners of irreducible representations of 
the octahedral group P0| . As the Vj(r) must be even against inversion, the Qj must have 



the same property, an only the six Qj with j = 1,6 discussed at the beginning of this section 
would appear in Eq. (fj), but we shall not consider the identical representation A\ because 
it does not modify the g tensor. The useful Vj(r) are given in Eqs. (4-6) from reference [JTO] . 



To study the effect that the V(r) given in Eq.(0) has on the g tensor of, we shall employ 
second order perturbation theory p5| , using both V(r) and the Zeeman term Hz = (g e S + 
L).H as perturbation. The change 5g in the g tensor is then obtained from 

S.5g.H = 2 - ( 9e + a)^ {-C E [V3Q 2 (S X H X - S y H y ) + 

Q 3 (3S £ H Z - S.H)] + C T [Qa (S z H y + S y H z ) + 

Qh (S X H Z + S Z H X ) + Q e (S x H y + S y H x )}} , (3) 

where \i& is the Bohr magneton, A is the splitting between the T 6 doublet and the lowest 
Tg quadruplet in the octahedral symmetry, and A = 283 cm" 1 in the P compound. 



The values of Ce and Ct are obtained by the same procedure employed in reference [10 



obtaining (r 4 ) from the cubic field parameter D q and (r 2 ) from the ratio (r 2 ) / yj (r 4 ) = 0.6544 

of the calculated values ^3|. For the Co-0 distance R we used R — 2.11176 A, corresponding 
to the reference octahedron defined below, and we found the values Ce = 6436 cm -1 /A and 
Ct = —3666 cm" - 1 / A. We can now calculate the crystal fields that would describe the 
experimental values of g or, what is equivalent, the corresponding normal modes within 
the approximations just discussed. As there are more normal modes than data, we fix the 
relations Q4 = Q3 = 0.3044Q6, which correspond to the normal modes calculated below 
from the crystallographic positions, and we obtain a perfect fit to the experimental values 
employing the normal modes given table [Vj. 

^From table p] we obtain the coefficients a = —0.9824, b = 0.1867, and a = —1.4128 
for Co2(OH)As0 4 . The ESR data was then adjusted with the normal modes coordinates 
given in table 0, where we used R = 2.1224 A, Q 4 = Q 3 = 0.4037Q 6 , C E = 6287 cm" 1 / A, 
Ct = —3558 cm -1 / A and A = 282 cm" 1 , for this compound. 

To calculate the crystallographic normal modes of the octahedron it is necessary to chose 
a reference perfect octahedron centered in the Co 2+ . To this purpose we consider the three 
normal modes of pure rotation, |Tj| Q 19 , Q 2 o, and Q21, and we chose the axes of the reference 
octahedron so that these three normal coordinates are zero, because they should not have 
any effect on the properties of the complex. The value R = 2.11176 A of the Co-0 distance 
in the reference octahedron was chosen so that Q\ = 0, and by this whole procedure we 
obtain a unique reference octahedron and minimize the effect of irrelevant rotations and 
expansions on the values of the normal modes. The direction cosines of the three Co-0 
directions in the reference octahedron are given in table [VI], and the normal modes derived 
from the crystallographic ionic positions given in table are shown in the third line of 
table |yj. The normal modes calculated from the crystallographic position of the O in the 
octahedral complex are different than those obtained from the experimental g tensor, given 
in the first line of the same table. This result indicates that although the nearest O to the Co 
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are the main source of the cubic field, Jl9| the remaining non-cubic perturbations have strong 
contributions due to the rest of the crystal. We conclude that the experimental g tensor 
could be explained by the crystal field V(r) of Eq. (Q) given in the axes of the reference 
octahedron defined in table [V]] with the Qj given in the row a) of table [V[ The agreement is 
perfect because there are more free normal coordinates than available 5g components, but 
the theory presented can only be considered a first approximation. In particular, although 
the crystal field theory of point charges gives the right symmetry properties, it is only a very 
rough description of the physics of the problem. 

Although we have not analyzed the hyperfine tensor in detail, we have verified that its 
components are compatible with the normal modes necessary to describe the 5g tensor. 

In the present calculation we have neglected the effect of the 4 T 2 triplet, that contributes 
to 5g in third order perturbation (our calculation would be of the second order). This effect 
was calculated by Tucker [23[] who obtained contributions that are about 6% of the second 
order contribution for the T 2 deformation and about 13% for the E deformation, and would 
therefore not alter substantially our conclusions. 



B. Penta-coordinated Co 



1 . The crystal field of the trigonal bipyramid. 



The structure of the penta-coordinated complex of Co 2+ is very close to a trigonal bipyra- 
mid, and the positions of the Co and of the five O are given with respect to the crystal axes 
in table [VI I|. Following a method similar to that employed in the octahedral case we chose 



an orthogonal system of axes X,Y,Z, such that the normal modes obtained from the crys- 
tallographic positions would not have contributions of irrelevant rotations and expansions. 
The direction cosines of the axes of this system with respect to the crystal axes a, b and c 
are given in table |V11I|, and the coordinates of the six atoms in the reference perfect trigonal 



bipyramid are given in table [TX|. There are two different Co-0 distances in the reference 
complex: R a corresponds to the three ligands in the XY plane (equatorial O) and R c to the 
two along the Z axis (axial O); their values for the phosphate and arsenate are given in the 
caption of table [TX| Two crystal field parameters D s and D t are necessary in the trigonal 
bipyramid, and are given in the point charge model by: ]lj ,2(| 



D s 



e 

14 

e 

168 



4g c 3g a 



(r 2 ), 



16g e 9% 

R 5 , R 5 „ 



(A 



(4) 



where we shall use q a = q c = ■ 
usual formula 

K/W = £ 

kq 



-2e. The crystal field potential V c t can be expressed by the 



47T 



2k + 



1 » r t> 



(5) 



where Yk q (6e, (pi) are the spherical harmonics at the position of the £-th ligand and the 
C^\6,ip) = \jAm I (2k + l)Y kq (9, (p) are usually called the Racah's rationalized spherical 
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harmonics. In our actual calculation we have employed the real combinations Ci m (8, ip) and 
Si m (0, p) that are proportional to cos(m, p) and sen(m, ip) respectively W^ - In the absence 
of the spin-orbit interactions one employs the irreducible representations V of the trigonal 
bipyramid to classify the eigenstates \a, S, L, V, 7, a) of the Hamiltonian, which are simply 
related to the states \a, S, L, Ml) (the index a identifies the particular states with the same 
S, L). In table I of reference |lTJ we find that the irreducible representations A 2 , A", A 2 , E' 
and E" are contained in the two terms 4 F and 4 P , and that the \a, S, L, Ml) states that gen- 
erate the corresponding subspaces are {|3, 3/2, 3, 0) , |3, 3/2, 1, 0)}— > A' 2 , { |3, 3/2, 3, ±3) } — > 
(A", A' 2 '), {|3, 3/2, 3, ±2)}^ E' and {|3, 3/2, 3, ±1) , |3, 3/2, 1, ±l)}-> E" . The Hamiltonian 
without spin orbit interaction is diagonal in the partners 7 of each irreducible represen- 
tation T and in the spin component Ms, so it is not necessary to write them explicitly 
here. The only C^\0,p) that contribute to Eq.(|5]) in the perfect trigonal bipyramid have 
k = 0,2,4 and q = 0. To calculate the matrix elements of the Hamiltonian that contains 
Vcf = 7 Vcf( r i), we have used the standard tensorial operator techniques as well as 
the unitary operators obtained form Nielsen and Koster's tables p8| , |29|| , and we have verified 
that our matrix coincides with that given in table II or reference [[Ll|] . 

Our main objective here is to find the gyromagnetic factors that one would expect to mea- 
sure in the penta-coordinated Co 2+ , and we shall employ the spectroscopic data measured 
by diffuse reflectance to estimate the parameters B, D s and D t for both Co:Zn 2 (OH)P0 4 
and Co:Mg 2 (OH)As0 4 . In the two rows labeled a) of table [H] we give the corresponding 
assignments of the transitions from the ground 4 A 2 to the levels with symmetry 4 A", 4 A 2 , 
4 E" , 4 E', 4 A' 2 (P) and 4 E"(P), where we use (P) to indicate the higher levels of the same 
symmetry. 

^From the eigenvalues of the Hamiltonian in the absence of the spin-orbit interaction, we 
find by trial and error the values of B, D s and D t that minimize the mean square deviation 
X for the two systems, and we give them in row b) of table [X[ The transitions calculated 
with these two sets of values are given in the two rows of table [TI] that are labeled b). The 
fitting is rather poor, and in particular the transitions to the levels 4 A", 4 A 2 and 4 E" fall 
below the range of the measuring equipment. As an alternative we have fitted only the three 
highest transitions, obtaining the values given in row c) of table [X|, and the corresponding 
values calculated with these two sets of parameters are given in the two rows of table that 
are labeled c). In the following section we shall consider these two sets of values to estimate 
the gyromagnetic factors for each of the two compounds. 



2. The spin- orbit interaction in the trigonal bipyramid. 

It is now essential to include the spin orbit interaction into the calculation. The ba- 
sis of the irreducible representations T 7 , T 8 and T 9 , of the double group D% h have a sim- 
ple expression in our system: they are given by \d 7 aSLJMj), and in particular we 
have r 7 (a) = {\d 7 aSLJ ± 1/2)}, r 7 (6) = {\d 7 aSLJ ± 11/2)}, T 8 (a) = {\d 7 aSLJ ± 5/2)}, 
r 8 (6) = {\d 7 aSLJ ±7/2)}, T 9 (a) = {\d 7 aSLJ ± 3/2)} and T 9 (b) = (\d 7 aSLJ ± 9/2)}. 
These states are easily obtained from the \d"' ', a, S, Ms, L, Ml) calculated above by employ- 
ing the 3-j or the Clebsch Gordan coefficients. In the absence of magnetic fields the two 
states of each Kramer's doublet have the same energy, and to calculate the energies of the 
system it is enough to consider only the states with positive Mj. As only the mixture of 



S 



the 4 F and 4 P states is important in our problem we shall consider only that subspace, 
and the corresponding matrix of the total Hamiltonian splits into five boxes of the following 
dimensions: (Mj = 1/2) -> T 7 (a) -»• (7 X 7), (Mj = 3/2) -> r 9 (o) ->(6x 6), (Mj = 
5/2) - r 8 (o) ^ (4 x 4), (Mj = 7/2) - r 8 (6) - (2 x 2), (M, = 9/2) - r 9 (6) - (1 x 1), 
and there are no matrix elements of Mj = 11/2, i.e. Tj(b), within the subspace { 4 F, 4 P} of 
d 7 that corresponds to S — 3/2. The matrices we have obtained coincide with those given 
in table II of reference [fL2|J , and their eigenvalues have been calculated for the different sets 
of B, D s and D t values that were obtained above, employing the one-electron spin-orbit pa- 
rameter C, = 580 cm -1 . For all the set of parameters in table |X| the lowest doublet is a T 7 (a) 
(Mj = ±1/2), separated by at least 75 cm -1 from the following T 9 (a) (Mj = 3/2) doublet, 
and by more than 2377 cm -1 from the remaining doublets. This situation is not altered by 
making fairly large changes in the three basic parameters B, D s and D t , and shows that 
even for moderate increases in the temperature only the lowest doublet (Mj = ±1/2) would 
be occupied. This doublet has allowed ESR transitions, and should be observed within the 
approximation employed. If the position of the two lowest doublets were exchanged, the 
ESR transitions of the lowest doublet would be forbidden and the spectra should not be 
then observed. 

The fact that the two lowest doublets have Mj = ±1/2 and Mj = ±3/2 and are separated 
by a large energy from the remaining doublets is easily understood when we notice that the 
lowest level in the absence of spin-orbit interaction is A A 2 . The orbital part A 2 is a singlet 
with no orbital angular momentum, and the total J would then correspond to the S = 3/2. 
These four states would be rather far apart from the remaining ones, and would split in the 
way calculated above through the higher order spin orbit mixing with those excited states. 

The present calculation was for a perfect trigonal bipyramid with D 3 h symmetry, and 
one wonders whether the deformations with respect to this structure could alter the relative 
position of the two lowest doublets, thus changing from an allowed to a forbidden ESR 
transition. We shall then study the effect of these deformations, both on the relative position 
of the two lowest doublets and on the value of the gyromagnetic tensor. In this study we shall 
follow a treatment similar to that employed in the octahedral case, by considering the effect 
of the normal modes of the trigonal bipyramid on the Hamiltonian of the penta-coordinated 
Co 2+ . 



3. The normal modes of the trigonal bipyramid 

As in the octahedral case we are interested in a contribution to the Hamiltonian of 
the same type of Eq. (|2|), but here the normal modes Qj and Vj(r) transform like the 
same partners of irreducible representations of the trigonal bipyramid. As the undistorted 
complex does not have a center of symmetry, both the even an odd modes against reflection 
in the equatorial plane may have non-zero matrix elements inside the configuration d 7 of 
Co 2+ , and therefore we shall need to consider both types of normal modes in our discussion. 

The departures of the six atoms of the complex span a reducible representation F of the 
D 3h group, that can be reduced as follows: T = 2A\ + A' 2 + 4E" ± 3A 2 + 2E" (see e.g. the 
Eq. (9.19) in reference |5D|). Of these irreducible representations, the A 2 corresponds to an 



axial rotation, one E" to two equatorial rotations, one A 2 to an axial translation and one 
E' to two equatorial translations. After eliminating these three translations and rotations 
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we are left with three even irreducible representations E', as well as two and one E" odd 
representations. The six even normal modes (Q\, Qe) transform in pairs like the partners 
of E'; they have been obtained employing standard techniques |30| , |3T| and are defined in 
table |Xl[ In the same way the two modes A' 2 ' (Qj, Q&) and the two partners of E" ((Qg, Q10) 
have been obtained, and are defined in the same table. 

The modes Q 3 , Q4, and Q 8 are translations of only the two axial oxigens, and the Q^.Qq, 
and Qi are translations of only the three equatorial oxigens while Qg and Qiq are rotations 
along the x and y axis of the two axial oxigens. As these modes are only partial rotations or 
translations they are capable of changing the crystal field. The x and y rotations of the three 
equatorial oxigens can be combined with Qg and Q10 to give full rotations of the trigonal 
bipyramid, and the z rotation of the three equatorial oxigens is already a full rotation, so 
these three sets of displacements would not appear in our calculation. 

/From tables |Vn|, |VlT3, and |X| we can calculate the displacements Uj of the five O 
with respect to their positions in the reference trigonal bipyramid (cf. section [Tl IB 1] ) and 
calculate the corresponding normal modes Qj defined in table |Xj. Employing these Qj we 
calculate in the next two sections the g values of the distorted complex. 



4- The effect of the normal modes on the crystal field 

We can now try and find an expression similar to Eq. (§) for the trigonal bipyramid. 
To this purpose we have employed a relation equivalent to Eq. (|5|) to calculate, for each 
of the ten normal modes Qj given in table [XI], the change in the crystal field V c f(r) when 
all the ligands are displaced from their equilibrium position in the reference complex by a 
small fraction e of that particular normal mode Qj. Expanding this change of V c f(r) in a 
power series of the coefficient e and taking the linear terms in e gives the corresponding 
V^ (r) from Eq. (||). As we are only interested in the subspace { 4 P, 4 P} with S = 3/2 of the 
configuration d 7 , and the Vj(r) are independent of the spin component M5, we need a 10x10 
matrix ( 4 L, M s , M L \V^ F \ 4 L', M s , M' L ) for each Q j: with fixed M s and L,L' = 3,1. There 
are regularities between the matrix elements associated to different Qj, and we shall employ 
the following abbreviations: Q a = (Q2 + 2Q1) +5(Q e + iQ 5 ), Qt = 3(Q 2 + iQi) + 7(Q 6 + iQ 5 ), 
Qc = 9(Q2 + iQi) + (Qe + iQs), Qd = iQi, and Q e = 9(Q W + iQg), as well as their complex 
conjugates Q*, Q* b Q* c , Q* d and Q* e . In tables [Xll| and [Xlll| we give the non-zero matrix 
elements in the upper triangle of the submatrices 4 F x 4 F and 4 P x 4 P respectively, and 
in table |X1 V| we give all those associated with 4 P x 4 P; the remaining non-zero elements 
are obtained by Hermitian conjugation. It is interesting to note that the matrix elements 
associated to Q 3 , Q 4 , and Q 8 are all zero: these modes involve only the two axial ions, and 
the corresponding two atom partial complex is not only invariant against the operations of 
D 3 h, but also against a twofold axis along the z direction. This extra symmetry forces all 
the one-electron matrix elements between d states of the crystal field associated to Q3, Q4, 
and Qs to be zero. 

As in the crystal field of the reference complex, the Vq F has coefficients containing the 
Co-0 distances R a and R c , as well as the atomic averages (r 2 ) and (r 4 ), and they appear 
as C2 and C4 in the tables |Xil| , |X1V| and |Xillf As the R a and R c are nearly the same, 



it is possible from Eq. (^) to relate the crystal field parameters D s and D t to these two 



10 



coefficients. Assuming that R a = R c we obtain c 2 = 14 D s and C4 = (168/25)D t , but we 
have derived slightly better relations considering the difference between R a and R c : 

7 

° 2 = -(3/2) + 2{R a /R c f Dsj 
21 

° A ~ (9/8) + 2 (R a /R c ) 5 ((>) 



As with the reference complex, we employ the 3-j coefficients to calculate the matrix 
elements of the crystal field V' CF in the representation that diagonalizes the total J and 
J z , because the doublets \d 7 aSLJ Mj) are basis for the irreducible representations of the 
reference trigonal bipyramid, and the eigenstates of the reference complex would then belong 
to subspaces with fixed Mj. In section [111 B 2| we have shown that, with the two sets of 
parameters B, D s and D t obtained in that section, the two lowest doublets belong to the 
Mj = 1/2 and Mj = 3/2 subspaces and that they are separated by more than 75 cm -1 , 
while the remaining doublets are more than 2300 cm -1 above them. A good approximation 
to calculate the effect of V'cp on these levels is then to consider the total Hamiltonian inside 
the two subspaces Mj = 1/2, 3/2, and one has then to consider a matrix of 26 x 26 elements, 
corresponding to values of J equal to 9/2, 1/2. The eigenstates of this matrix show that 
there is no change in the relative position of the two lowest doublets, so that the ground 
state remains Mj = 1/2. 



5. The g-factors of the penta- coordinated Co 2+ . 



To calculate the spin Hamiltonian we employ the traditional method |25fl . In the present 
case we consider the four states of the two lowest doublets of the reference trigonal bipyramid 
calculated in section [D 1B2| as the eigenstates of the unperturbed Hamiltonian, with Mj = 



1/2 as the ground doublet and Mj = 3/2 as the excited one. Both the Zeeman term and the 
crystal field V' CF produced by the deformation of the normal modes are the perturbations, 
and in the usual way we find the gyromagnetic tensor g in second order. We have calculated 
the three components of g for the penta-coordinated Co 2+ for all the sets of B, D s and D t 
given in table |X|, and the results are given in table X.V. The values corresponding to the 



reference trigonal bipyramid are given in the rows &o an d Co, while those given in the rows bj 
and Cj (with j = 1,2 for the phosphate and j = 1,2,3 for the arsenate) have been calculated 
employing the normal modes Qj derived from the crystallographic positions (cf. table |V11| ) 
as discussed in section II 1 B 3| . One verifies in table |XV| that the average of the principal 



values of g is not very different from 4.33, but that it changes with B and the crystal field 
parameters more than in the octahedral case. The rows b\ and c\ include the effect of all 
normal modes, while in b 2 and c 2 only the even modes are considered. In the phosphate 
case, the crystallographic Qi = Q§ = 0, and for the arsenate we have also imposed this 
condition in rows 63 and C3. /,From table |XV| we conclude that 

• With only even modes and Qi = Q5 = (compare rows b 2 and c 2 in phosphate and 
63 and C3 in arsenate with rows b\ and ci) the average g is not altered by the lower 
symmetry crystal fields generated by the remaining normal modes, and these fields 
affect the equatorial components of g, but leave their sum and the axial component 
unaffected. 
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• Only the axial g-factor is altered by the inclusion of the odd modes, while the two 
equatorial ^-factors are not altered (compare rows 61 with 6 2 and c% with C2). This 
result is true both with Qi = Q 5 = (in the phosphate) or otherwise (in the arsenate). 

• The fields associated to the modes Qi and Q$ change the sum of the two equatorial 
^-factors but leave the axial value unaffected (compare rows 62 with 63 and C2 with C3 
in the arsenate). 

^From the calculations in the present section, it follows that one should observe an 
allowed ESR line of Co 2+ from the penta-coordinated complex when that site is occupied. We 



have seen in Section |II Q that besides the lines associated to the octahedral spectra, there are 
some weak extra lines that could be interpreted as belonging to that complex: their estimated 
g-factors are given in rows b) of table [III], and should be compared with the values given 



in table XV, that were calculated for different sets of parameters derived from the optical 



spectra and with normal modes calculated from the crystallographic positions. It is clear that 
the arsenate values in row C3 of table pCV] are fairly close to the estimated values in row b from 



table |TJ. It is well known that the g(i) obtained from the crystallographically calculated 
normal modes are generally different from those experimentally observed, as discussed for 
the octahedral compounds (cf. section 111 A| ), and we could expect that a good fitting 



could be obtained by making small changes in the crystallographic normal modes. To verify 
this assumption it is sufficient to change only Q2 and Qq, keeping all the remaining modes 
at their crystallographic values. Employing Q2/R a = —0.03 and Q%jR a = —0.07 for the 
phosphate we find g(l) = 7.05, g{2) = 3.03, and g(3) = 2.12, while for the arsenate we 
obtain g{l) = 7.62, g(2) = 3.04, and g(3) = 1.99 with Q 2 /R a = 0.01 and Q 6 /R a = -0.07. 
These fittings are fairly good, and show that the ESR spectra of Co 2+ in the two compounds 
can be described perfectly well within our theory, but the crystal fields obtained can not be 
taken too seriously because of the very large errors in the experimental g tensor. 

We notice that the relative intensities of the extra lines in Fig. |2| are rather smaller than 
those in Fig. |^. This can be understood because the concentration of Co 2+ in the arsenate is 
ten times larger than in the phosphate, and this should also alter their relative occupations. 

The rather small intensity of the lines that could be attributed to the penta-coordinated 
complex indicates a very small occupation of Co 2+ in the penta-coordinated sites. To ver- 
ify this conjecture, we present a molecular orbital calculation of the heat of formation of 
these compounds in the following section, and the results are compatible with the present 
conclusion. 



C. Molecular Orbital Calculations 

The energetics of penta- and hexa-coordinated phosphate clusters has been assessed in 
terms of the molecular orbital theory. Heats of formation were calculated within the well 
known semi-empirical technique Parametric Model 3 (PM3) |32|]. This is a technique derived 
from the Hartree-Fock approximation in combination with a minimal basis set expansion 
of the molecular orbitals. Here we used a special parametrization developed for transition 
metal atoms which is contained in the package SPARTAN Correlation and relativistic 
effects, which are not explicitly treated in this theory, are partly recovered from the adoption 
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of experimental data in the parametrization. Molecular geometries were obtained as follows: 
the central metal atom and the coordinates of the first neighboring five or six oxygen atoms 
were taken from the crystal structure of the compound Co2(OH)P04. The ligands were 
chosen to be phosphoric acid molecules, OP(OH) 3 , since they have all bonds saturated and 
are neutral. Geometry optimizations of the isolated ligand were carried out at the ab-initio 
6-3 1G** level of calculation, in which the atomic orbitals of the basis set are written as a 
linear combination of cartesian gaussian functions . The proton-free oxygen atoms of the 
ligands were then placed in the crystallographic positions of the oxygens around the metal 
atom such that the 0=P bond points to the M-0 direction, as shown in Figs. |] and [5]. The 
PM3 heats of formation of the clusters [M(OP(OH) 3 ) n ] 2+ , n = 5,6 and M = Co and Zn were 
calculated assuming that Zn ion just replace the Co ion at frozen ligands positions. This is a 
reasonable assumption since the pure Co and Zn crystals have very similar cell parameters. 
Co clusters are doublets so that the unrestricted (spin-polarized) PM3 Hamiltonian was 
adopted. Spin contamination was negligible in this calculation. In order to discount the 
energies associated to the ligands themselves, the heats of formation of the corresponding 
clusters without the central metal ion were computed. Results are shown in table [XV1| . The 
values in the first column are the contributions from the ligands to the metal clusters heats 
of formation. It is then possible to evaluate the relative stabilization of Co 2+ and Zn 2+ 
ions in the penta- and hexa-coordinated environments by making the difference between the 
values in columns two or three and column one. It gives the energies —1258.04 kcal/mol and 
— 1314.30 kcal/mol, for Co 2+ at the trigonal bipyramid and octahedral sites respectively, 
while for Zn 2+ the values are 346.19 kcal/mol and 340.00 kcal/mol. These values show 
that Co prefers the octahedral site by an amount of ~56 kcal, which is approximately 2.4 
eV, and Zn is slightly more stable also at the octahedral site by ~6 kcal, or 0.3 eV. This 
difference is due to the partially filled 3d orbitals of Co that interact with the lone pairs of 
the neighboring oxygens, giving a more covalent character to the interaction. 

A more direct comparison was made through the calculation of the heat of forma- 
tion of clusters where the metal atoms and first neighbors are in the conformation of the 
Co 2 (OH)P0 4 unit cell, as illustrated in Fig || Two clusters were built such that in one the 
Co 2+ ion occupies the hexa-coordinated site and Zn 2+ ion is in the penta-coordinated site 
(cluster 1), while in the other (cluster 2) Co 2+ and Zn 2+ ions are exchanged. The phosphate 
ions in contact with both metal ions were substituted by H2PO4 species and the remain- 
der of the ligands were phosphoric acid molecules The composition of these clusters is then 
[CoZn(OH)(H 2 P0 4 )2(OP(OH)3) 5 ] 4+ . PM3 calculations gave AH f (cluster 1) = -2147.10 
kcal/mol and AHf (cluster 2) = —2094.44 kcal/mol, that is, the cluster with Co 2+ ion in the 
octahedral site is ~53 kcal more stable than the other one. It is thus expected that Co impu- 
rities in the zinc compounds occupy preferentially the octahedral sites, and this conclusion 
agrees with very low intensity of the ESR lines attributed to Co 2+ in the penta-coordinated 
sites of the dilute compounds, as discussed in the previous section. 

IV. DISCUSSION AND CONCLUSIONS 

Four compounds of the adamite type: Zn 2 (OH)P04, Mg 2 (OH)As0 4 , Co2(OH)P04, and 
Co 2 (OH)As0 4 .have been synthesized and studied, and the measurement of the optical prop- 
erties of the pure Co compounds and of the ESR of impurities of Co 2+ in Zn 2 (OH)P04 and 
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Mg 2 (OH)As0 4 have been discussed. Crystal field theory has been employed to try and un- 
derstand the experimental ESR results for the two Co 2+ complexes with coordination five 
and six that are present in the adamite structure. The Racah parameter B as well as the 
crystal fields D q for the octahedral complex and both D s and D t for the trigonal bipyramid 
one have been estimated from the assignments that were made of the diffuse reflectance 
spectrum of these two complexes. Two alternative sets of parameters were proposed for the 
penta-coordinated complex. 

iFrom the crystallographic structure, a reference octahedron centered in the Co 2+ was 
defined, such that the normal modes of the complex corresponding to rotations and expan- 
sions would be zero and the remaining normal modes would not have any contribution of 
these irrelevant deformations. Using a method already applied [pi]] to study the ESR of 
Co 2+ in NH4MPO4.6H2O, the crystal fields that would reproduce the experimental g tensor 
of the octahedral complex in both Zn 2 (OH)P04 and Mg 2 (OH)As0 4 have been obtained. 

As the penta-coordinated complex seems to have at most minor contributions to the 
ESR spectra, we have analyzed the possible motives for this behavior. We argue that two 
doublets with Mj = ±1/2 and Mj = ±3/2 would be lowest in energy, separated by a rather 
large excitation energy from the remaining excited states. The Mj = ±3/2 has forbidden 
ESR transitions, and this would explain the experimental results if that were the ground 
doublet, but when the crystal field of the trigonal bipyramid is considered together with 
the spin-orbit interaction, it was found that the Mj = ±1/2 doublet, with allowed ESR 
transitions, is the lowest. To verify whether this result would be altered by the deformations 
of the trigonal bipyramid, we considered their effects in a way similar to that employed in 
the octahedral case to calculate the g tensor. First it was necessary to derive the normal 
modes of the trigonal bipyramid that are relevant to our problem, and they are given in table 
XI . The corresponding Jahn- Teller contributions V' CF to the crystal field, whose non-zero 
matrix elements ( 4 L, M s , M L \ V' CF \^L' , M s , M' L ),foi L, L' = 3, 1, are given in tables [XTl], |XTV 
and |X.111| . Defining a reference perfect trigonal bipyramid by the same method employed in 
the octahedral case, the values of the relevant normal modes were obtained by employing 
the crystallographic positions, and subsequently used to calculate their effect on the relative 
position of the two lowest doublets. No appreciable change was found, and as an alternative 
explanation we assumed that the penta-coordinated complex is scarcely occupied in the 
dilute system. To verify this conclusion, the heat of formation of the octahedral and the 
trigonal bipyramid complexes with both Co and Zn as the central ions were calculated, 
and it was found that their values are compatible with a rather small occupation of the 
penta-coordinated site. 

Employing the Jahn- Teller crystal fields together with the normal modes calculated from 
the crystallographic distortions, it was possible to calculate the g tensor, shown in table 
XV| for both the perfect and deformed trigonal bipyramid, this last subjected to different 
deformations. The trace of the g tensor in the perfect trigonal bipyramid changes more with 
the parameters B,D S and D t than in the octahedral case, where it is always fairly close to 
13. 

The trigonal bipyramid has no center of symmetry, and it was necessary to consider all 
the normal modes, even those that are odd against reflection in the horizontal symmetry 
plane. We have shown that these last modes affect the axial component of g but that 
they have little or no effect on the two equatorial components. The experimental g-tensor 
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of the penta-coordinated complex could be measured but with rather large errors. As in 
the octahedral case, the crystallographically determined normal modes could not explain 
the observed values, but for the two type of complexes it was possible to find values of the 
normal modes that would generate crystal fields that describe the experimental ESR spectra 
for both the phosphate and arsenate compounds. 

We conclude that both our theoretical analysis of the ESR of the systems studied, as well 
as the molecular orbital calculation of the formation energies, coincide in assigning a rather 
small relative occupation of the penta-coordinated sites with respect to the octahedral ones 
in those systems. We could also understand the experimental spectra of both the octahedral 
and penta-coordinated complexes by considering the effect of the crystal fields generated by 
the corresponding normal modes on the g tensor. 
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TABLES 





4 T 2g 


4 A 2s 


4 Ti 9 (P) 


B 




P0 4 


a) 


8450 


15450 


18350 






b) 


7819 


16013 


18324 


767.6 


819.4 


As0 4 


a) 


7700 


15500 


18020 






b) 


7616 


15585 


18011 


758.9 


796.9 



TABLE I. a) The transitions between the ground 4 Ti 9 level and the levels shown at the top 
of each column, in cm -1 and assigned from the experimental spectra of the octahedral complexes 
of Co2(OH)P0 4 and Co2(OH)As0 4 ; the level A T\ g {P) corresponds to the highest of the same 
symmetry, b) The best fit, obtained with the B and D q shown in the last two columns. 



Ji 2 — > 


4 A" 4 A" 

^H) ^2 


4 E" 


4 E' 


*A' 2 (P) 


4 E"(P) 


P0 4 


a) 
b) 
c) 


6400 
3233 
1511 


7000 
4835 
3603 


11100 
12868 
11106 


15800 
17386 
15801 


19600 
17947 
19604 


As0 4 


a) 
b) 

c) 


5000 
2707 
1437 


6250 
4440 
3535 


10870 
12210 
10876 


16000 
17188 
15999 


19800 
18595 
19805 



TABLE II. a) The transitions between the ground level 4 A" and the levels shown at the top of 
each column, given in cm -1 and assigned from the experimental spectra of the penta-coordinated 
complexes of Co2(OH)P0 4 and Co2(OH)As0 4 . b) The best possible fit to the five transitions, c) 
The best fit obtained by adjusting only the three transitions of higher energy. The corresponding 
values of B, D s , Dt are given in rows b) and c) of table [X|. 



gi g2 g3 Ai A 2 A 3 









P0 4 








a 


5.89±0.02 


4.55±0.05 


2.02±0.02 


240±5 


155±8 


85±3 


b 


8. ±0.5 


3.2 ±0.3 


2.0 ±0.2 














As0 4 








a 


6.22±0.02 


4.21±0.05 


2.05±0.02 


140±5 


120±7 


55±5 


b 


9. ±1.5 


3. ±0.5 


2.0 ±0.2 









TABLE III. Values of the principal g and A parameters, obtained from the spectra in Figs. |2j 
and m The values of the A parameters are in 10 -4 cm -1 units, a) Octahedral complex: the g and 
A values were obtained from a program simulating powder spectra, as described in the text, b) 
The parameters for the penta-coordinated Co 2+ , also estimated by simulation. 
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11 cx 1 _M. ^ 




C (^AJ 


a (a; 




C (A) 










ASO4 




1 4.7150 


-1.2420 


0.0000 


0.9197 


-1.0794 


0.0000 


2 4 9627 


— 1 0327 


2 9700 

— ■ >) \ \)\) 


6438 


—1 2633 


3 0255 


3 2.1592 


-1.2713 


1.7274 


-1.9091 


-1.1662 


1.3476 


4 3.0793 


1.0327 


2.9700 


-0.6438 


1.2633 


3.0255 


5 3.3270 


1.2420 


0.0000 


-0.9197 


1.0794 


0.0000 


6 5.8827 


1.2713 


1.7274 


1.9091 


1.1662 


1.3476 


7 4.0210 


0.0000 


1.5135 


0.0000 


0.0000 


1.4910 


TABLE IV. The columns a,b,c give 


the position of the six oxyj 


;ens (n=l,...,6) and Cobalt 


(n=7) in the hexa-coordinated complexes 


of Co 2 (OH)P0 4 


and Co2(OH)As0 4 with respect to the 


three unit cell axes. The X,Y,Z axes roughly correspond to 


n=l,2,3 respectively, taking 


n=7 as the 



origin. 



Q 2 /R Q 3 /R Qa/R Q5/R Qe/R 



P0 4 


a) 





-0.04407 


0.03478 


0.03478 


0.11425 


b) 





-0.11683 


-0.01717 


-0.01717 


-0.05641 


As0 4 


a) 





-0.0441 


0.04032 


0.04032 


0.09989 


b) 





-0.09328 


-0.02818 


-0.02818 


-0.06982 



TABLE V. Normal modes of the octahedral Co divided by the Co-0 distance R in 
Co:Zn2(OH)P0 4 and Co:Mg2(OH)As0 4 . a) Values that adjust the experimental values of the 
g tensor, b) Values obtained from the crystallographic positions corresponding to the pure com- 
pounds. 



X Y Z X Y Z 







P0 4 






As0 4 




a 


0.4054 


-0.5794 


-0.7071 


0.3800 


-0.5963 


-0.7071 


b 


0.4054 


-0.5794 


0.7071 


0.3800 


-0.5963 


0.7071 


c 


-0.8194 


-0.5733 


0.0000 


-0.8434 


-0.5373 


0.0000 



TABLE VI. Direction cosines of the three axis X,Y,Z of the reference perfect octahedron of 
Co2(OH)P0 4 and Co2(OH)As0 4 with respect to the three crystallographic axis a, b, c 
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n 


a (A) 
a [A) 


h (A\ 
a (A) 


r (A) 
C [A) 


a (A) 
a {A) 


h (A) 


r (A) 
C (A) 






i ^ / \ 

PO4 






A / \ 

ASO4 






4 9008 


3 2438 


2 9700 

_■-./( \J\J 


4 9744 


91 32 


3 01 80 


2 


1.8617 


2.9132 


1.2426 


1.9003 


1.1604 


1.3442 


3 


1.8617 


2.9132 


4.6974 


1.9003 


1.1604 


4.6676 


4 


3.1412 


5.1252 


2.9700 


3.2085 


3.2015 


3.0180 


5 


3.0793 


1.0327 


2.9700 


3.2736 


-0.9132 


3.0180 


6 


2.9102 


3.0621 


2.9700 


3.0056 


1.1561 


3.0180 




TABLE VII. The columns a,b,c give 


the position of the five oxyj 


$ens (n=l,.. 


,5) and Cobalt 


(m 


=6) in the penta- 


coordinated complexes 


of Co 2 (OH)PO, 


t and Co2(OH)As0 4 with respect to the 



three unit cell axes. 



X Y Z X Y Z 







P0 4 






As0 4 




a 


0. 


0. 


1. 


-0.00203 


0. 


0.999998 


b 


0.99931 


0.03720 


0. 


0.999698 


-0.02449 


0.00203 


c 


-0.03720 


0.99931 


0. 


0.02449 


0.9997 


0.00005 



TABLE VIII. Direction cosines of the three axis X,Y,Z employed to define the reference perfect 
trigonal bipyramid of Co2(OH)P0 4 and Co2(OH)As0 4 with respect to the three crystallographic 
axis a, b, c 



n X Y Z 

1 

2 —^R a —\Ra 

3 -^-R a ^ 

4 R c 

5 -R c 

6 

TABLE IX. The columns X,Y,Z give the position of the five oxygens (n=l,...,5) and Cobalt 
(n=6) in the reference perfect penta-coordinated complexes of Co2(OH)P0 4 [R a = 2.01622 A and 
R c = 2.04365 A], and Co 2 (OH)As0 4 [R a = 1.98578 A and R c = 2.05596 A] with respect to the 
axes defined in table VQI| 
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B 


D s 


D t 


B 


D s 


D t 






P0 4 






As0 4 




b 


728. 


165. 


947. 


785. 


313 


919. 


c 


852. 


745. 


885. 


875. 


749. 


869. 



TABLE X. The values of B, D s , Dt in cm" 1 that fit the optical transitions, given in table ||, 
of the two penta-coordinated complexes. The best fit to the five transitions is given in row b), and 
the best fit to the three highest transitions in row c). The spin orbit parameter £ = 580 cm -1 was 
used in all these fittings 





Xl 


yi 


Zl 


x 2 


Y2 


Z2 


X3 


y3 


Z3 


x 4 


y4 


z 4 


X5 


ys 


Z5 


2\/3Qi 


-2 








1 


V3 





1 


V3 























2V3Q2 





2 





V3 


-1 





V3 


-1 























V2Q3 





























1 








1 








V2Q4 
































1 








1 





V3Q 5 


1 








1 








1 


























VZQs 





1 








1 








1 























V3Q 7 








1 








1 








1 




















V2Q 8 



































1 








1 


V2Q 9 
































-1 








1 





V^Qio 





























1 








-1 









TABLE XL The six even (1-6) and four odd (7-10) normal modes E' that are relevant to our 
problem. The numbers are the coefficients of the departures Uj = {xj,yj, Zj} of the j=th ion from 
their equilibrium position. 
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Mr 

±VJ - Li 


Mr 

L 


( 4 _F Mq Mr WAtp^F M? Mr >) 


—6 


o 

— Z 


— CZ + 1U C4J Q e f [i\OJ 


g 

— 6 


— 1 


(0 C2 Q a — 15 C4 g?, J / (ODV 5 J 




n 




-3 


1 


-\/5 c4 Q*/16 


-2 


-1 


-(9 c2 - 40 c4) Q*/ (21\/5") 


-2 





(12 c2 Q a + 5 c4 Q 6 ) /(56>/T0) 


-2 


1 


v/5 c4 Q rf / C 2a/6^ 


-2 


2 


-5 c4 Q*/(16\/3) 


-1 





- ^2 (3 c2 - 25 c4) Q*/ (35\/3~) 


_i 

i- 


1 


r9 O -1-95 r4 Di.^ /('14D\/3 N I 


.\ 


2 


-a/5 r4 Oj/ ( 2\/fi^ 


_1 


3 


-\/5 r4 O* /lfi 





1 


a/2 (3 c2 - 25 c4) Q*/ f 35\/3) 





2 


(12 c2 g a + 5 c4 Q 6 ) / (56>/T0) 





3 


-^15 c4 Q d /4 


1 


2 


(9 c2 - 40 c4) Q*/ (21^/5) 


1 


3 


(6 c2 g a - 15 c4 Q b ) /(56\/5) 


2 


3 


(3 c2 + 10 c4) Q*/ (7y/3) 



TABLE XII. The non-zero matrix elements of the crystal field V' CF = J2j Qj ^j'( r ) generated 
by the normal modes Qx, ...,Q\q between states \a, S, Ms, L, M£) = | FMsMl) in the subspace 
4 i ? x 4 F. Only the elements corresponding to the upper triangle of the matrix are given, and the 
remaining ones are obtained by Hermitian conjugation. The matrix is independent of, and diagonal 
in, the spin components Mg- To compress the table we have used the following abbreviations: 

Q a = (g 2 + 1Q1) + 5(g 6 + iQ 5 ), Qb = 3(g 2 + iQi) + 7(g 6 + *g 5 ), g c = 9(g 2 + %Qi) + (g 6 + iQ 5 ), 

(even modes) and Qd = 1Q7, Q e = (1Q9 + Q10) (odd modes), as well as their complex conjugates 
Ql Qt Qt, Q* d and Qt 



M L 


M' L 






CP, M s , MlIV^P, M s , M' l >) 


-1 
-1 





1 
1 




3 c2 g*/5 
-3y/3 c2 Q a /20 
-3 c2 g*/5 


TABLE XIII. 


Same as in table 


XII 


but for the sub-matrix 


4 P x 4 P. The same abbreviations 



are used here. 
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M L 


M' L 


< 4 P,M S ,M L |V£ F | 4 F,M S ,M£) 


-1 


-3 


(72 c2 Q* -5 c4 Q£)/(56^30) 


-1 


-2 


\/2(12c2 + 5 c4) Q e /(7y/T&\ 


-1 





-2(18 c2 + 25 c4) Q*/105 


-1 


1 


(24 c2 Q a - 25 c4 Q 6 ) / (280^) 


-1 


2 


c4 Q d /4 


-1 


3 


-V5 c4 QH (8x/6) 





-3 


-V§ c4 Q d / Uy/t) 





-2 


(12 c2 Ql + 5 c4 Q*) / (28>/l0) 





-1 


\/2 (24 c2 - 25 c4) Q e / (35\Z3j 





1 


-y/2 (24 c2 - 25 c4) Q e / (35\/3) 





2 


(12 c2 Q a + 5 c4 Q 6 ) / (28^Io) 





3 


-y/E c4 Q d / (4^3) 


1 


-3 


-V5 c4 Q c / (8>/6) 


1 


-2 


c4 Q d /4 


1 


-1 


(24 c2 Ql - 25 c4 Q*) / (280^2) 


1 





2 (18 c2 + 25 c4) Q*/105 


1 


2 


-V2(12 c2 + 5 c4) Q e / (7\/l5) 


1 


3 


(72 c2 Q a - 5 c4 Q 6 ) / (56\/3o) 



TABLE XIV. Same as in table XD but for the sub-matrix 4 P x 4 F. All the non zero matrix 
elements are given here, and those corresponding to the sub-matrix 4 F x 4 P are obtained by 
Hermitian conjugation. The same abbreviations are used here. 
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gl 


g2 


g3 


Sav 


P0 4 




4.8027 


4.8027 


1.9904 


3.8653 


h 


5.0477 


4.5577 


2.2998 


3.9684 


b 2 


5.0477 


4.5577 


1.9904 


3.8653 


CO 


5.0435 


5.0435 


1.9829 


4.0233 


Cl 


5.7379 


4.3492 


2.1171 


4.0680 


C2 


5.7379 


4.3492 


1.9829 


4.0233 


As0 4 


60 


4.8723 


4.8723 


1.9885 


3.9110 


h 


5.7637 


4.2279 


2.2364 


4.0760 


b 2 


5.7637 


4.2279 


1.9885 


3.9934 


b 3 


5.6402 


4.1044 


1.9885 


3.9110 


CO 


5.0667 


5.0667 


1.9818 


4.0384 


Cl 


7.1232 


3.5385 


1.9854 


4.2157 


C2 


7.1232 


3.5385 


1.9818 


4.2145 


C3 


6.8591 


3.2744 


1.9818 


4.0384 



TABLE XV. The principal components of the calculated g tensor for the penta-coordinated 
Co 2+ and their average g av in Co:Zn2(OH)P0 4 and Co:Mg2(OH)As0 4 . The rows bo and cq are for 
the reference trigonal bipyramid, while bi, b 2 , &3, ci, C2, and C3 include the effect of deformations 
produced by the crystallographically calculated normal modes. Rows b\ and c\ include all the 
normal modes, while 62 an d c 2 only include the even modes, and in rows 63 and c e we have also 
put Q\ = Q 2 = 0. The values of B, D s , D t and C employed here for rows bj ( j = — 3) are given 
in row b of table [X], and those corresponding to rows Cj are given in row c of that table. 



AH j (kcal/mol) 


no ion 


Co 2 + 


Zn 2 + 


[M(OP(OH) 3 ) 5 r 


-1214.69 


-2472.73 


-868.50 


[M(OP(OH) 3 ) 6 ] 2+ 


-1440.16 


-2754.45 


-1100.15 



TABLE XVI. Heats of formation, in kcal/mol, from PM3 calculations. 
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FIGURES 



FIG. 1. Diffuse reflectance spectra for the a) Co2(OH)P04 and b) Co2(OH)As0 4 . The horizon- 
tal scale is linear in the wavelength, but has been labeled employing the corresponding wavenum- 
bers. 

FIG. 2. ESR spectra of Co 2+ in Zn2(OH)P04. a) Experimental spectrum, b) Sum of the 
simulated spectrum for both the hexa-coordinated and penta-coordinated complexes, c) Simu- 
lated spectrum for the hexa-coordinated complex. The g-values of both the octahedral and of 



the penta-coordinated complexes are given in table ID. The arrows in c) show the g-values and 
their positions for the octahedral complex, while those in a) correspond to the penta-coordinated 
complex. The insert gives the detail of the experimental hyperfine structure and of the simulated 
one (around 52)) attributed to the Co 2+ in the triangular bipyramid. 

FIG. 3. ESR spectra of Co 2+ in Mg2(OH)As0 4 . Curves a), b), and c), and the meaning of the 
arrows and values in curves a) and c) are the same as in figure ^. 

FIG. 4. Ball-and-stick model of the penta-coordinated metal clusters with phosphoric acid 
molecules as ligands. One atom of each type is labeled in the figure. 

FIG. 5. Hexa-coordinated metal/phosphoric acid cluster. One atom of each type is labeled in 
the figure. 

FIG. 6. Cobalt/Zinc cluster built from the atomic coordinates of the Co2(OH)P04 unit cell. 
Part of the phosphate ions were replaced by phosphoric acid molecules. One atom of each type is 
labeled in the figure, which shows the Co in the penta-coordinated position 
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